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ABSTRACT 

We construct axisymmetric Schwarzschild models to measure the mass profile of the local group dwarf 
galaxy Fornax. These models require no assumptions to be made about the orbital anisotropy of the stars, as 
is the case for commonly used Jeans models. We test a variety of parameterizations of dark matter density 
profiles and find cored models with uniform density pc = (1-6 ±0.1) x 10"^ M© pc"^ fit significantly better than 
the cuspy halos predicted by cold dark matter simulations. We also construct models with an intermediate- 
mass black hole, but are unable to make a detection. We place a 1-cr upper limit on the mass of a potential 
intermediate-mass black hole atM.< 3.2 x 10'^ Mq. 

Subject headings: dark matter — galaxies: dwarf — galaxies: individual (Fornax) — galaxies: kinematics and 
dynamics — Local Group 
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1. INTRODUCTION 

Low-mass galaxies provide a unique testing ground for pre- 
dictions of the cold dark matter (CDM) paradigm for struc- 
ture formation, since they generally have a lower fraction of 
baryons than massive galaxies. These galaxies allow for a 
more direct measurement of the underlying dark matter po- 
tential, as the complicated effects of baryons on the dark 
matter are less pronounced. A particularly testable predic- 
tion of CDM is that all galaxies share a universal dark mat- 
ter density profile, characterized by a cuspy inner power law 
p oc r"" where a = I (Navarro et al. 1996b, hereafter NFW). 
Many authors have investigated low-mass spirals and found, 
in contrast to the predictions of CDM, dark matter density 
profiles with a fla t inner core of slope a = ( Burkertjjl995t 
Persic et al." 19961; 'de Bl oket all 120011: iBlaTs-OueUette et all 
2001; Simon et al. 2005). This has launched the debate 
known as the core/cusp controversy. 

A number of other studies have investigated the mass con- 
tent o f dwarf spheroidal galaxies (dSphs). iGilmore et al.l 
(l2007h give a comprehensive review of recent attempts to con- 
strain the i nner slope of their dark matter profiles with Jeans 
modeling ('Jeans"1919'; 'Binne y & Tremainal 19871 chapter 4). 
When significant, cored profiles are preferred for all dSphs 
modeled (Gilm ore et alJI2007l and references therein). 

These results, however, are subject to a major caveat of 
Jeans modeling; it is complicated by the effect of stellar ve- 
locity anisotropy. Models fit to the line-of-sight component of 
the velocity dispersion, but anisotropy can severely affect the 
modeling of enclosed mass. Therefore, additional assump- 
tions must be made. The studies presented in Gil more et al, 
(2007) assume spherical symmetry and isotropy. Evans et al I 
(12009) show that a weakness of Jeans modeling is that given 
these assumptions combined with the cored light profiles ob- 
served in dSphs, the Jeans equations do not allow solutions 
w ith anything o ther than a cored dark matter profile. 

'Wa lker et al.l ('2009b) construct more sophisticated models 
and attempt to parameterize and fit for the anisotropy. As a re- 
sult, preference for cored profiles becomes model-dependent. 
They therefore are unable to put significant constraints on the 
slope of the dark matter profile. This highlights the main prob- 
lem with Jeans modeUng — it is highly dependent on the as- 



sumptions made. 

Distribution function models are more general than Jeans 
models, and progress has been made applying them to a 



ap pr 

number of dSph^ systems (iKlevna et al.f B002r IWul l2007t 



lAmorisco & Evansll201 Ibl) . Nevertheless these models still 
make strong assumptions such as spherical symmetry or 
isotropy, and models that do fit for anisotropy do so without 
using the information about the stellar orbits contained in the 
line-of-sight velocity distributions (LOSVDs). 

We employ a fundamentally different modeling technique, 
known as Schwarzschild modeling, that allows us to use this 
information to self-consistently calculate both the enclosed 
mass and orbital anisotropy. Schwarzschild modeling is a 
mature industry, but one that has seldom been applied to the 
study of dSph galaxies (see Valluri et al. 2005(). 

In addition to being well-suited for measuring dark mat- 
ter profiles, Schwarzschild modeling has often been used to 
search for black holes at the centers of galaxies. Another un- 
resolved issue relevant to the study of dSphs is whether they 
host an intermediate-mass black hole (IMBH). In a hierar- 
chical merging scenario, smaller galaxies are thought to be 
the building blocks of larger galaxies. It is thought that all 
massive galaxies host a supermassive black hole (SMBH) at 
their center, therefore it is logical to believe that their build- 
ing blocks host smaller IMBHs. Evidence for these IMBHs is 
scarce, however, and dynamical detections are even scarcer. 
The closest and lowest mass example of a dynamical mea- 
surement is an upper limit on the local group dSph NGC 205 
ofM.<2.2x IO^Mq dValluri et al.l2005!) . Black holes in this 
mass range can provide constraints on theories of black hole 
growth and formation. The two most prominent competing 
theories of nuclear black hole formation are direct collapse 
of primordial gas (jU memura et al. 1993t lEisenstein & Loebl 
ll995tlBegelman etal. 2006) or accretion onto and mergers of 
seed black holes re sulting from the collapse of the first stars 
dVolonteri & Pernall2005l) . 

In this paper we present axisymmetric, three-integral 
Schwarzschild models in an effort to determine the inner slope 
of the dark matter density profile as well as the orbit structure 
of the Fornax dSph. We also investigate the possibility of a 
central IMBH . We assume a distance of 135 kpc to Fornax 
(lBersieril2000l) . 
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Fig. 1 . — Line-of-sight velocity distributions of four bins. Open circles with 
error bars are the data. Over-plotted are the model values for the best-fitting 
cored model (red) and NFW model (blue). Bins are located at: (a) R = 297" , 
e= 18° (b)R = 550",6»= 18° (c) R = 1008", 6» = 45° (d) R = 2484", 5 = 45°. 
Quoted values are un-reduced. 

2. DATA 

To construct dynamical models, we require a stellar light 
profile as well as stellar kinematics in the form of LOSVDs. 
We use published data for both the photometry and kinemat- 
ics, and describe the steps taken to convert this data into useful 
input for our models. 

2.1. Stellar Density 

To determ ine the stellar density, we use a number density 
profile from lColeman et alJ (120051) extending to 4590". We 
linearly extrapolate the profile out to 6000" — a physical ra- 
dius of 3.9 kpc at our assumed distance. We also extrapolate 
the profile inwards at constant density from 90" to 1". 

To convert to a more familiar surface brightness profile we 
apply an arbitrary zero-point shift in log space, adjusting this 
number so that the integrated pr ofile ret urns a luminosity con- 
sistent with the value listed in iMateol (!l998). Adopting an 
ellipticity of e = 0.3 (lMateol ll998'), we deproject under the 
assumptio n that surfaces of con stant luminosity are coaxial 
spheroids (iGebhardt et al.|[T996h . and for an assumed inclina- 
tion of ! = 90°. 

2.2. Stellar Kinematics 

We de rive LOSVDs from in dividual stellar velocities pub- 
lished in I Walker et al.l (l2009al) . The data contain heliocentric 
radial velocities and uncertainties with a membership proba- 
bility for 2,633 Fornax stars. Most of these are single-epoch 
observations, however some are multi-epoch. Stars that have 
more than one observation are averaged, weighted by their 
uncertainties. After making a cut in membership probabil- 
ity at 90%, we are left with 2,244 stars. Although a signif- 
icant number of stars observed may be in binary or multiple 
systems, simulations have shown that such systems are un- 
likely to affect measured dispersions tHargreaves et al...l996t 
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Fig. 2. — Gauss-Hermite moments for stars near the major axis (blue), mi- 
nor axis (red), and averaged over all angles (green). Solid lines correspond 
to the best-fit model with a cored dark matter halo, dashed lines are for the 
best-fit model with a NFW halo. 



[Olszewski et al.lll996tlMateolll998l) . ^ 

We adopt a position angle PA = 41° ( Wal ker et all 12006*). 
We assume symmetry with respect to both the major and mi- 
nor axes and fold the data along each axis. To preserve any 
possible rotation, we switch the sign of the velocity whenever 
a star is flipped about the minor axis. 

The transverse motion of Fornax contributes a non- 
negligible line-of-sight velocity to stars, particularly those 
at large galac tocentric radius. U sing the equations in Ap- 
pendix A of IWalker et alJ (120081) . we correct for this ef- 
fect. We adopt values for t he proper motion of (/Xq.,A*5) = 
(47.6,-36.0) mas century"' (" Piatek et al.l 120071) and assume 
the heliocentric r adial velocity of Fornax is 53.3 km s~' 
dPiatek et al.ll2002h . 

We divide our meridional grid into 20 radial bins, equally 
spaced in approximately log r from l"to 5000". There are 5 
angular bins spaced equal ly in sin 6* over 90° from the major 
to the minor axis (G ebhardt etalj|2000t ISiopiset al] 12009) . 
From the positions of the folded stellar velocity data, we de- 
termine the best binning scheme so that each grid cell contains 
at least 25 stars from which to recover the LOSVD. Our first 
bin with enough stars to meet this criterion is centered at 47", 
and the last bin is centered at 2500". We therefore have two- 
dimensional kinematics coverage over the radial range 47"- 
2500"(30pc - 1.6 kpc). At small radii the number density 
of stars with velocity measurements is low, thus our central 
LOSVDs have higher uncertainty compared to those at larger 
radii. 



Within each grid cell, we calculate the LOSVD from dis- 
crete stellar velocities by using a n adap tive kernel density 
estimate adapt e d from ISilvermanI ( 119861) and explained in 
IGebhardt et al.l ( 119961) . We estimate the 1 - cr uncertainties 
in the LOSVDs through bootstrap resamplings of the data 
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TABLE 1 
Best-Fit Model Parameters 



DM Profile 5: c r, (kpc) pdMQpc-^) M.(Mq) N,„m 

NFW 239.8 1.3±0.6 4.1±0.26 I1.7±1.4 — — 3124 

Log 162.6 1.5±0.5 — — 1.6±0.1 x 10"^ — 4319 

Log 162.6 1.6 ±0.2 — — 1.6±0.1 x 10"^ < 3.2 x lO'* 3423 



Note. — Best-fit parameters for NFW, and cored logarithmic dark matter halos. is un-reduced, the 
number of degrees of freedom are the same for each model. Model parameters and l-cr uncertainties are 
quoted. Wmorfe/ lists the number of models run for the corresponding parameterization. 



(IGebhardt et al.l[T996HGebhardt & Thomasll2009h . The boot- 
strap generates a new sample from the data itself by randomly 
picking data points, where is the total number of stars in 
a given bin, allowing the same point to be chosen more than 
once. We then estimate the LOSVD from that realization and 
repeat the procedure 300 times. The 68% confidence band 
on the LOSVDs corresponds to the 68% range of the realiza- 
tions. We compare the velocity dispersion as measured by the 
LOSVDs with the biweig ht scale (i.e., a ro bust estimate of 
the standard deviation, see Beers et al.|[T990l) of the individual 
velocities and note good agreement. 

Figure [U plots the LOSVDs of four bins. Rather than pa- 
rameterizing these LOSVDs with Gauss-Hermite moments, 
our models instead fit directly to the LOSVDs to constrain the 
kinematics of the galaxy. However, we do fit Gauss-Hermite 
moments for plotting purposes only. These data are presented 
in Figure |2] for stars that have been grouped into bins near 
the major axis (blue) and minor axis (red). Near the center 
of the galaxy the density of stars with kinematics is sparse, 
so we therefore group stars into annular bins covering all an- 
gles (green). We estimate the I-ct uncertainties of the Gauss- 
Hermite moments by fitting to each of the 300 realizations 
calculated during the bootstrap discussed above. The error 
bars plotted contain 68% of the 300 realizations. 

3. DYNAMICAL MODELS 

The modeling code we use is d e scribe d in detail in 
pebhardt et al. (2003).ThomaseLalJ (12004 pb05), and 
[Siopis et al. (2009) and is based on the technique of orbit 
superposition ( Schwarzschild 1979). Similar axisym metric 
codes ar e described in Rix et al. ( 1997); van der Mar el et alj 
1X998"); 'Cretto n et ah 0999^; .Valluri et al., (.2004) while 
Ivan den Bosch et al.l (120081) present a fully triaxial 
Schwarzschild code. Our code begins by choosing a 
trial potential that is a combination of the stellar density, 
dark matter density, and possibly a central black hole. We 
then launch ^ 15,000 orbits carefully chosen to uniformly 
sample the isolating integrals of motion. In an axisymmetric 
potential, orbits are restricted by three isolating integrals of 
motion, E, Ly, and the non-classical "third integral" /?. As 
it is not possible to calculate h a priori, we use a carefully 
designed sch eme to systematically sample h for e ach pair 
of E and L. (iThomas et al.ll2004t ISiopis et al.ll2009h . Orbits 
are integrated for many dynamical times, and each orbit is 
given a weight w,. We find the combination of w, that best 
reproduces the observed LOSVDs and light profile via a 
minimization subj ect to the constraint of maximum entropy 
(ISiopis et al.ll200si) . 

We run models by varying 3 parameters — the stellar MjLy 
and two parameters specifying the dark matter density profile. 
Some models are also run with a central black hole whose 
mass is varied in addition to the other 3 model parameters. 



Each model is assigned a value of and we identify the best- 
fitting model as that with the lowest x^. We determine the 
68% confidence range on parameters by identifying the por- 
tion of their marginalized x^ curves that lie within Ax^ = 1 of 
the overall minimum. 

3.1. Model Assumptions 

Our trial potential is determined by solving Poisson's equa- 
tion for an assumed trial density distribution. On our two- 
dimensional polar grid, this takes the form: 



M 

p{r, 0) = —^{r, 0) + pDMir) 



(1) 



where M/L is the stellar mass-to-light ratio, assumed constant 
with radius, and v(r, 6) is the unprojected luminosity density. 
The assumed dark matter profile poMif) is discussed below. 
For simplicity, we assume Fornax is edge-on in all our mod- 
els. 

3.2. Dark Matter Density Profiles 

We parameterize the dark matter halo density with a num- 
ber of spherical density profiles. We use NFW halos: 



PDuir) = 



200 A(c)p„,, 
3 (r/r,)(l + r/r,)2 



(2) 



where 



A(C): 



ln(l + c)-c/(l+c) 



and Peril is the present critical density for a closed universe. 
The two parameters we fit for are the concentration c and scale 
radius r,. We also use halos derived from the logarithmic po- 
tential: 



PDM(r) ■■ 



3rg + r^ 
4^(r2 + r2)2 



(3) 



These models feature a flat central core of density pc = 
yv} /AiiGrl for r <rc and an r~^ profile for r > re- We fit 
for Vc and Tc, the asymptotic circular speed at r = cx) and core 
radius respectively. We run over 10,000 models with only 
three distinct parameterizations: NFW halos, and logarithmic 
models with and without an IMBH. 

4. RESULTS 

We find significant evidence for cored logarithmic dark 
matter density profiles. These models are preferred at the 
Ax^ = 77 level when compared to models with an NFW 
halo, a highly significant result. Perhaps more convincingly, 
the values for the concentration preferred by our models are 
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Fig. 3. — x' curves for all parameterizations of the mass profile. NFW halos (blue) ai'e parameterized by concentration c and scale radius r,. Logaiithmic halos 
with an IMBH (green) and without (red) are specified by Vc and r^. We also plot core density pc = 3V'c^/47rGr^ as it is the controlling parameter over the radial 
range of our models. We fit for stellar MjLy in all models (upper left panel). NFW models have much higher and are scaled down by 75 to fit on the same 
axis. Black hole mass for logarithmic halos with an IMBH (green) is plotted in the upper right panel. Note the apparent minimum in for logarithmic halos with 
an IMBH is due to incomplete parameter sampling. 



around c = 4. Only relatively recently formed structures like 
galaxy clusters are expected to have concentrations this low 
dNFW,) . 

Table[T|summarizes the results of our models, while Figures 
[T|and|2]illustrate the preference for cored models over models 
with an NFW halo in fitting to the kinematics. We stress again 
that LOSVDs like those plotted in Figure [T] are the kinematic 
constraint, and not the Gauss-Hermite moments of Figure|2] 

While we fit for Vc and in the cored models, these pa- 
rameters are strongly degenerate. Our model grid extends 
to 3.3 kpc, thus any model with > 3.3 kpc has a uni- 
form density pc = 3V^ /AirGr^ over the entire range of our 
model. Furthermore, we have no velocity information from 
stars past R>1.6 kpc and therefore cannot constrain the kine- 
matics in the outer parts of the galaxy. Thus, for models with 

> 1 .6 kpc, Pc is now the only parameter that differentiates 
between models. As pc is dependent on both V^. and r^, the 
latter two parameters are completely degenerate. 

Figure |3] illustrates this effect. Plotted are the curves 
for each model parameter. Lines of the same color indicate 
a common parameterization of the mass profile (e.g. cored -i- 
IMBH). While the for both Vc and asymptotes to large 
values, Pc is tightly constrained. Note that the behavior of 
rc for logarithmic profiles with an IMBH (green line) is a re- 
sult of incomplete parameter sampling. With a more densely- 
sampled parameter space, the curve for rc for cored models 
with an IMBH would likely asymptote to large r^. in a similar 
fashion as models without an IMBH (red curve). 

The addition of a central black hole to the mass profile does 



not make a noticeable difference to the overall x^ for most 
values of M,. We therefore place a 1-cr upper limit on M, < 
3.2 X IO'^Mq. 

We plot the mass profile for our best-fit model in Figure 
|4](solid black line with surrounding 68% confidence region). 
This is a cored logarithmic dark matter profile without a cen- 
tral black hole. The mass profile of our best-fit dark halo is 
plotted as the dashed line, and the stellar mass profile is plot- 
ted in red. The contribution of dark matter to the total mass 
increases with radius as the local dynamical mass-to-light ra- 
tio rises from approximately ^ 2 to greater than 100 in the 
outermost bin of our model. 

4.1. Orbit Structure 

We construct a distribution function for the galaxy from the 
set of orbital weights w, resulting from the x^ minimization 
of our best-fit model. To explore the orbit structure, we deter- 
mine the internal (unprojected) moments of the distribution 
function in spherical coordinates. Streaming motions in the 
r and 6 directions are assumed to be zero. In this coordinate 
system, cross-terms of the velocity dispersion tensor are zero. 

Figure |5] plots the anisotropy in the diagonal components 
of the dispersion tensor While some panels show an average 
value near unity, there are regions in every panel where the 
ratio plotted is different from one. Additionally, we define 

the tangential velocity dispersion a, = ^ \^{{v]f,) + cr^) where 
(v^) is the second moment (v^) = (t^ + V^, and is the mean 
rotation velocity. With this definition, we plot the ratio ar/a, 
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Fig. 4. — Total enclosed mass for our best-fit model (black line with sur- 
rounding confidence region). Red line is the enclosed stellar mass. Dashed 
line is our best-fit dark matter halo. 

in the bottom panels of Figure|5]to investigate whether orbits 
are radially or tangentially biased. From these plots it is clear 
that the common assumptions of Jeans modeling — constant or 
zero anisotropy — are unrealistic. We find that at most radii in 
the galaxy, orbits are radially biased. The uncertainty in the 
anisotropy is largest at small radii, as evidenced by the size 
of the 68% confidence regions in Figure |5] This is likely due 
to the sparsity of kinematics in the inner part of the galaxy 
(there are limits to how closely target fibers can be spaced in 

multi-fiber spectroscopy)^ 

In a recent paper, Kaza ntzidis et al.l (1201 1[) simulated the ef- 
fects of tidal stirring on a number of dSph progenitors around 
a Milky Way sized halo. They found radial anisotropy in all 
of the final remnants, and our models are consistent with these 
findings. 

5. DISCUSSION 

5.1. Cores and Cusps 

Our analysis shows that for the Fornax dwarf an NFW dark 
matter halo with inner slope a = 1 is rejected with high confi- 
dence. We have kinematics from 30 pc-1 .6 kpc, and over this 
range the models prefer an a = uniform density core with 

= 1.6 X \Qr^ Mq^c~^ . We do not attempt to fit for models 
with an intermediate value of the slope < a < 1 . Further in- 
vestigation is necessary before we can conclude that the best 
fitting dark matter profile is the logarithmic model. The steep 
a = 1 cusp of the NFW profile is, however, robustly ruled out. 

The models, in general, seem to prefer less mass in the areas 
over which we have kinematic constraints. In NFW models, 
the concentration c sets the normalization (or y-intercept) of 
the density profile. Because c cannot be lowered below an 
astrophysically reasonable limit, NFW models enclose more 
mass than cored models. This difference is reflected in the 
difference between cored and NFW models, as the kinematics 
are best fit by models with less mass. Figure |2] hints at this 
as the best fit NFW model (dashed line) typically has higher 



values for a than either the data or best-fitting cored model 
(solid line). 

Several groups have approached the core/cusp issue in 
dSphs by taking advantage of the fact that some dSphs 
host multiple populations of tracer stars that are chemi- 
cally and dynamically distinct. By fitting models to each 
component, the underl ying dark matter profile can be mod- 
eled more accurately. lAmorisco & EvansI (1201 lal) fit two- 
component distribution fun ction models to Sculptor, while 
IWaUcer & Penarrubial (1201 1 *) apply a convenient mass estima- 
tor (discussed below) to each stellar component in Sculptor 
and Fornax. It is believed that this mass estimator is unaf- 
fected by orbital anisotropy, thus their method yields a robust 
determination of the dynamical mass at two locations in the 
galaxy — allowing for the slope of the dark matter profile to 
be measured. Each of these studies finds models with a cored 
dark matter halo preferable to the predicted cuspy NFW pro- 
file. 

It must be noted, however, that we are not observing the 
pristine initial dark matter distribution in this galaxy. Rather, 
it has likely been modified by complex baryonic processes 
over the lifetime of the galaxy. These processes may in- 
clude: adiabatic compression (Blumenthal et al. 1986), halo 
rebounding fo llowing baryonic mass loss from supernovae 
dNavarro et al.l [T996a) . or possibly dynamical friction act- 
ing on clumps of baryo ns (lEl-Zant et al.ll200Tt but see also 
iJardel & Sellwoodll2009l) . Although we chose this galaxy be- 
cause these effects were likely to be small, they are neverthe- 
less not well understood and our result must be taken in that 
context. 

5.2. Central IMBH 

We are unable to place a significant constraint on the mass 
of a central IMBH. Figure|3](upper right) shows the marginal- 
ized curve against IMBH mass for cored dark matter den- 
sity profiles. The curve asymptotes to low values of IMBH, 
thus we are only capable of placing an upper limit on the mass 
of any potential IMBH. Furthermore, our best-fit cored model 
with and without an IMBH have the same y^. We therefore 
impose a 1-ct upper limit on M, < 3.2 x 10*Mq. It is unfor- 
tunate that we are not able to place a lower limit on M, be- 
cause measurements of black holes in the range M. <104 Mq 
place direct constraints on SMBH formation mechanisms (van 
Wassenhove et al. 2010 ). Our models, however, do robustly 
rule out a black hole of larger mass. 

In massive galaxies it is thought that the radius of influence, 
^inf ^ GM, I must be resolved in order to detect and pre- 
cisely measure a black hole (G ebhardt et al. 2003; Kormend^ 
l2004t iFerrarese & FordI l2005rTGultekin et al.il2009h . Using 
our upper limit on M, we can calculate the maximum radius 
of influence of a potential black hole. Estimating the central 
velocity dispersion at ct ~ 10 km s~' gives an upper limit for 
^inf < 14 pc. Our kinematics start at = 26 pc, so it is not sur- 
prising that the minimum black hole mass we were able to de- 
tect has Rinf close to 26 pc. To detect smaller black holes, we 
require kinematics of stars closer to the center of the galaxy. 

We are able to detect the dynamical influence of a black 
hole with a similar mass as Valluri et al. (2005) detect in 
NGC 205, however with kinematics of much lower resolution. 
Our innermost model bin is centered around 30 pc whereas 
they use high-resolution kinematics from the Hubble Space 
Telescope to resolve spatial scales less than 1 pc. The advan- 
tage we have is that the central velocity dispersion is much 



6 



1.5 
1.0 
0.5 
0.0 
1.5 
1.0 
0.5 
0.0 
1.5 
1.0 
0.5 
0.0 
1.5 
1.0 
0.5 






100 1000 100 1000 

r (arcsec) 



Fig. 5. — Anisotropy in various components of the velocity dispersion ten- 
sor. Sliaded regions correspond to tlie 68% confidence regions, solid lines 
plot the best fit model. Left and right hand panels plot stars near the major 
and minor axes, respectively. 



smaller in Fornax, which makes Ri^f larger for fixed M,. 
NGC 205 is also more than five times as distant as Fornax. 

5.3. Mass Estimators 

Several authors have come up with convenient estimators 
of total mass within a given radius for local group dSphs. 
[Strigari et al. (2008) use the mas s enclosed wit hin 300 pc 
while Wal ker et al.l i2009b) and W olfetal.1 dloTO) find a sim- 
ilar expression for the mass contained within the projected 
and un-projected half-light radii, respectively. These es- 
timators bear striking resemblance to a result obtained by 
iCappellari et al.l ( 120061) derived from integral field kinematics 
of massive elliptical galaxies, and they all hint at an easy way 
to determine dynamical masses without expensive modeling. 
They are believed to be insensitive to ve locity anisotropy 
based on the derivation in lWolf et al.l (120101) . and we compare 
their estimates to our models as a check on this. 
For the mass contained within 300 pc we measure M300 = 
X lO^MfTi, roughly a factor of three smaller than 



o c+0.77 
■-'•-'-0.11 
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iStrigari et al] (120081) who measure Msoo = 1-14!^: 12 
parameterized anisotroi 
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using Jeans mo dels with pararn et erizea anisotr o py. 

The^ CappellimetalJ d2006l) . IWalker et all (l2009bb . and 
IWolf et al.i (12OIOI) mass estimators are all of the form: 



-0.1 1 



X 10'Mq and 
group's value for 
M(R,) « 
8.1!^i X 
lO'^M, 



o 



only by the value of the constant k. In order to more fairly 
compare between these estimators and our models, we use 
the values for the luminosity-weighted line-of-sight velocity 
dispersion {(jIos) ~ ll-3;!;[ gkm s~' projected half-light radius 
Re = 689pc, and un-projected half-light radius re = 900pc that 
we calculate from the data used in our models. 

Our best-fitting model has M{Re) = 3.Q+"46 
Mire) = 5.8!(^:*] X 10^ Mq. With each 
k and our kinem atics, the mass estim ates are 
5.l!i:g X IO^Mq (IWalker etaH l20()9bh . M(r,) 
10^ Mr ("Wolfet al.' '20l^, and M(i?,) « 1.0^ ^ 
(ICaDDellari et al. 2006). Our mode l is broadly consi stent with 
both the iWalker et all (I2009bl) and IWolf et al.1 dMTOI) estima- 
tors. 

The evidence that mass estimators are anisotropy- 
independent comes l argely from comparis on to spherical 
Jeans models (except ICappellari et al.ll2006l) . The weakness 
of these models is that the anisotropy must be parameterized 
and is restricted to be a function of radius only. Our models 
are not subject to these constraints since the anisotropy is cal- 
culated non-parametrically and is free to vary with position 
angle. We suggest that the best way to prove the accuracy 
of mass estimators is to compare with models that can self- 
consistently calculate both mass and anisotropy for realistic 
potentials. 

For bright elliptic al galaxies, ICappellari et al.l ( |2006l) and 
iThomas et al.l (1201 ll) have done just that. In these cases, the 
mass estimates are checked against masses derived from ax- 
isymmetric Schwarzschild modeling and good agreement is 
found. Ours is the first study to perform a similar test with 
dSphs, and there is no reason to assume that success with 
bright ellipticals guarantees accuracy in the dSph regime. The 
results from our comparison above are nevertheless reassur- 
ing. 

5.4. Tidal Effects 

The principle of orbit superposition, and hence our entire 
modeling procedure, relies on the assumption that the galaxy 
is bound and in a steady state. The amount of tidal stripping in 
Fornax due to the effect of its orbit through the Milky Way's 
halo is not well-known. For reasonable values of Fornax total 
mass m. Milky Way mass M, and Galactocentric radius /?o, the 
tidal radius of Fornax is r, ~ {m/3M)^I^RQ ^13.5 kpc. This 
estimate of r, is sufficiently larger than our model grid that we 
would not expect tidal effects to be important over the radial 
range of our models. If Fornax is on an eccentric orbit about 
the Milky Way, however, the above equation for r, is not valid 
and estimation of the tidal radius is not as straightforward. 
Fortunately, studies investigating its trans verse motion sug- 
gest the orbit of Fo rnax is roughly circular (iPiatek et alJl2007l : 
IWalker et al.ll2008h . 



M(re J = k{alos)Re 



(4) 



where rest is the radius at w hich the estimator is v alid. For 
ICappellari et all ( l2006l) and IWalker et all (l2009bl) = Re 
(the projected half-light radius), while for lWolf etall (120101) 
'"est = re (the un-projected half-light radius). Other than 
the projected/un-projected difference, each estimator differes 
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